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Abstract. The Strong Equivalence Principle (SEP) de- 
mands, besides the validity of the Einstein Equivalence 
Principle, that all self-gravitating bodies feel the same ac- 
celeration in an external gravitational field. It has been 
found that metric theories of gravity other that than gen- 
eral relativity typically predict a violation of the SEP. In 
case of the Earth-Moon system (weak field system) this 
violation is called the Nordtvedt effect. 

It has been shown by Damour and Schafer, that 
small-eccentricity long-orbital-period binary pulsars with 
a white dwarf companion provide excellent conditions to 
test the SEP in strong field regimes. 

Based on newly discovered binary pulsars this paper 
investigates a possible violation of the SEP in strong field 
regimes. New limits with an improved confidence level are 
presented. 

The results of this paper lead to constrains on the com- 
bination s/2+( of the only two (post)^-Newtonian param- 
eters s and C that arise from the (post) ^-Newtonian ap- 
proximation of the tensor-multi-scalar theory of Damour 
and Esposito-Farese. 
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1. Introduction 

When dealing with relativistic theories of gravity one is 
confronted with three types of equivalence principles (see 
Will 1993): 

— the Weak Equivalence Principle (WEP), 

— the Einstein Equivalence Principle (EEP), 

— and the Strong Equivalence Principle (SEP). 

The WEP, which goes back to Galileo and Newton, states 
the universality of free fall for all neutral test masses. In a 
more geometrical view the WEP states that all test masses 



move along geodesies in spacetime. Test masses are under- 
stood to be bodies with negligible self-energy and there- 
fore with negligible contribution to spacetime curvature. 
The WEP has been verified with a fractional precision of 
better than 10~^^ (Roll et al. 1964, Braginsky & Panov 
1971, Su et al. 1994). 

The EEP demands, besides the validity of the WEP, 
that in local Lorentz frames the non-gravitational laws of 
physics are those of special relativity. The EEP implies 
that spacetime has to be curved and thus is the basic 
ingredient of any metric theory of gravity. 

The SEP states, besides the validity of the EEP, the 
"universality of free fall for self-gravitating bodies" . One 
has to be careful with the notion of a freely falling self- 
gravitating bodies in an external gravitational field. There 
is no rigorous definition for the SEP in relativistic theories 
of gravity. Because of non-linearity the split of the metric 
field into an external and a local part can only be approx- 
imate. For a discussion of the SEP within a slow-motion 
weak-field approximation see Ashby & Bertotti (1984) and 
Bertotti & Grishchuk (1990). For metric theories of grav- 
ity, other than general relativity, it has been found that 
they typically introduce auxiliary gravitational fields (e.g. 
scalar fields) and thus predict a violation of the SEP (see 
e.g. Will 1992, Will 1993). 

In 1968 Nordtvedt (1968a,1968b) proposed to test a 
possible violation of the SEP through the analysis of lunar 
laser ranging (LLR) data. In case of a violation of the SEP, 
the external field of the sun would cause a "polarization" 
of the lunar orbit since the Earth and the Moon show 
different fractions of self-energy. A violation of the SEP 
implies the that the ratio of passive gravitational mass 
nig and inertial mass m,- differs from one by a function of 
the gravitational self-energy ii/srav q£ ^j^g body: 



nii 



1 + A[E8rav] 
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Nordtvedt parameter. Present LLR data lead to 



-0.0016 <T]< 0.0006 
-0.0008 <r]< 0.0013 



(Dickey et al. 1994), 
(Miiller et al. 1995). 



(2) 



Damour and Schafer (1991) pointed out that, in view 
of the smallness of the self-gravity of planetary bodies 
(e.g. £'8''™/mic2 = -4.6 x 10-^° for the Earth), such 
solar-system tests of the SEP indicate nothing about 
higher-order gravitational-energy contributions to the ra- 
tio rrig/mi. 

For binary pulsars the situation is different, since for a 
neutron star one finds E^™ /rriiP ~ —0.2, which takes us 
into the strong field regime. Especially small-eccentricity 
long-orbital-period binary pulsars with a white-dwarf 
companion [E^™ /rriiP ~ -10"'*) constitute an excellent 
laboratory to test higher-order effects in the violation of 
the SEP. 

Only recently, Damour and Esposito-Farese (1995) 
presented a new, field-theory-based framework for dis- 
cussing and interpreting experimental tests of relativis- 
tic gravity, notably at the (post)^-Newtonian order. They 
use a model in which gravity is mediated by a tensor field 
together with one or several scalar fields ( "tensor- multi- 
scalar" theories, see Damour & Esposito-Farese 1992). 
Damour and Esposito-Farese come to the conclusion that 
within their framework the (post)^-Newtonian deviations 
from general relativity can be fully described by intro- 
ducing two new parameters, e and C, beyond the usual 
(Eddington) parameters /3 = /3 — 1 and 7 = 7 — 1 of the 
standard parameterized post-Newtonian (PPN) formal- 
ism (see Will 1993). They further find that the parame- 
ters e and C are practically immeasurable by present and 
near future solar-system experiments. They conclude that, 
because of the importance of self-gravity effects in neu- 
tron stars, binary-pulsar experiments are an ideal testing 
ground for the (post) ^-Newtonian structure of relativistic 
gravity, e.g. small-eccentricity long-orbital-period binary 
pulsars with a white-dwarf companion can be used to set 
limits on the combination e/2 -|- C. 

In this paper we investigate newly discovered binary 
pulsars (Taylor & al. 1993, Taylor et al. 1995) with regard 
to a violation of the SEP. In Sect. 2 we display the theoret- 
ical background for testing the SEP in small-eccentricity 
long-orbital-period binary pulsars. Sect. 3 introduces the 
test systems. In Sect. 4 we present the results and in Sect. 
5 we give the conclusions. 

2. Violation of the SEP and binary motion 

In this section we repeat briefiy the calculations given by 
Damour & Schafer (1991). The masses of pulsar and com- 
panion will be denoted by rrip and rric, respectively. The 
total mass of the binary system is given by M = rrip + rric. 
The ratio of passive gravitational mass and inertial mass 
for pulsar and companion will be denoted by 1 -|- Ap and 



The orbital period of the binary system will be denoted 
by Ph and the orbital eccentricity by e. The eccentricity 
vector e is a vector directed along the apsidal line towards 
the periastron. 

In case of a violation of the SEP the equations of rel- 
ative motion r{t) have the form 



QM- 



aR + Ag, 



(3) 



where Q denotes the effective gravitational constant for the 
interaction between rrip and rric, slr is the orbital relativis- 
tic corrections of post-Newtonian order, and g denotes the 
gravitational acceleration field of the Galaxy, g is practi- 
cally constant in magnitude and direction, so that we can 
interpret A g as a "gravitational Stark effect" caused by a 
possible violation of the SEP.** 

For small-eccentricity long-orbital-period binary pul- 
sars Eq. (3) implies that the orbital plane is fixed and 
that the eccentricity vector e{t) is given by the following 
vectorial superposition: 



e{t) = GA -I- eR{t). 



(4) 



The constant vector ga is directed along the projection of 
g onto the orbital plane (gi). 
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ujR is the average angular velocity of the relativistic ad- 
vance of the periastron, n = 27r/_P5, and a is the semi- 
major axis of the relative motion, ga represents a con- 
stant, g-induced, polarization of the orbit. eR{t) is a vec- 
tor of constant length which lies in the orbital plane and 
rotates with angular velocity ujr, and it represents the 
usual relativistic advance of the periastron. Figure 1 illus- 
trates the temporal behaviour of the eccentricity vector 




Fig. 1. Evolution of the eccentricity vector e(t) = Ga -|- e_R(t) 



It is obvious that the limit of validity of Eq. (4) is that ujr 
should be appreciably larger than the angular velocity of 



** We assume that the limits for a violation of the WEP ob- 
tained on Earth apply also for white-dwarf and neutron-star 
matter. 



reference frame of the binary system. 

Allen & Martos (1986) developed a mathematically 
simple Galactic mass model that gives a good fit to the 
observed Galactic accelerations. This model will be used 
here to determine g at any given point in the Galaxy. 
In the Galactic plane we find then for circular orbits the 
following approximation: 



O'.'04/yr 
p[kpc] 



for 



P ■ 



1 



lOkpc, 



(6) 



where p denotes the distance to the Galactic center. Most 
of the binary pulsars considered below are rather close 
to the Galactic plane, i.e. Eq. (6) can be used safely to 
estimate the ratio ujr/uj for most of our test systems. For 
ujR one finds 



c'^a{l — e^) 



(7) 



where in general relativity the factor T is unity and Q = 
GNewton = G- For small-eccentricity binary pulsars in the 
solar vicinity (p « 8 kpc) one finds = I, Q = G): 



ujr/u) « 1.4 



Ph 



lO^days 



-5/3 / ^ ^ 2/3 



(8) 



Thus only for pulsars with -P6 ^ 10^ days we expect that 
we cannot use the simple solution Eq. (4). 

For a small-eccentricity binary pulsar with observed 
eccentricity e and ujr being appreciably larger than uj, 
the inequality 



l/sm9 for 9 e [0, Itt) 

L2-l^] (9) 



1 for 6* e \^TT, Itti 



- 1 / sin 6 for 



holds (see Fig. 1). Using the orbital inclination, i, the lon- 
gitude of the ascending node, SI, and the angle between 
the direction of sight and g. A, one obtains from Eq. (9) 
the following limit on the violation of the SEP: 



\A\<-f,,xie,n), 

e 



where 



and 



m 



[1 — (cos i cos A -|- sin i sin A sin Sl)^]!/^ 



(10) 
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the binary system one gets with rather high precision 
directly from pulsar timing. 

— g and A can be calculated using the Galactic- 
acceleration model of Allen & Martos (1986). These 
two quantities depend on the location of the binary 
pulsar in the Galaxy and thus on the direction of sight 
and on the distance Earth-pulsar. The distance of a 
pulsar is calculated from the observed dispersion in 
the signal, based on measurements of the density of 
the interstellar distribution of free electrons (Taylor 
& Cordes 1993). The error in distance is typically of 
the order of 25%. The error in the direction of sight is 
extremely small. 

— The inclination i and the total mass M can be deter- 
mined if one knows the mass of the pulsar nip and the 
mass of the companion rric: 

M = m, + m„ sm* = — (13) 

where x = Gp sini/c is the (observable) projected semi 
major-axis of the pulsar orbit. 

The analysis of neutron-star neutron-star binary pul- 
sars implies that with 95% confidence the lower 
bound for neutron-star masses is between I.OIMq and 
I.34M0 and the upper bound is between 1.43M0 and 
1.64M0 (Finn 1994). Since during the mass transfer 
a considerable amount of mass could have accreted 
onto the neutron star (Bhattacharya & van den Heuvel 
1991) the upper limit for our test systems could ex- 
ceed 1.64M0. The combined investigation of X-ray bi- 
naries and radio-binary pulsars (Thorsett et al. 1993, 
van Kerkwijk et al. 1995) suggests I.8M0 as an (con- 
servative) upper limit for the mass of a pulsar. Thus 
we will use I.OM0 < nip < I.SMq. 
The mass of the companion white dwarf is limited by 
theoretical studies of stellar evolution in binary stars 
(Joss et al. 1987, Savonije 1987, Rappaport et al. 1995). 
For small-eccentricity binary pulsars nic is given as a 
function of the orbital period Ph- 

— In case of the small-eccentricity binary pulsars PSR 
J1713-F0747 and PSR B1855-F09 (will not be used 
below) it was possible to measure the Shapiro delay 
caused by the companion. By this we get further lim- 
its on nip and nic (Camilo et al. 1994, Kaspi et al. 
1994). 

— The angles 6 and SI are not observable. We shall treat 
them as independent random variables, uniformly dis- 
tributed between and 27r. Thus, one can only give 
limits on |A| with a certain probability by excluding 
those areas of the parameter space where fi \ is 
larger than some given value. 

3. The test systems 



On the right-hand side of Eq. (10) we have eight quanti- 
ties, six of which can be deduced from observations. 



To be able to apply Eq. (10) we need binary pulsars which 
show the following properties: 



nary system with a large value for P'^ / e. As a selecting 
criterion for useful test systems one finds / e > 10^ 
days^. (Including binary systems with Pj^/e < 10^ 
days^ would not change the results of this paper, since 
they put only weak limits few %) on |A|.) 
— The binary pulsar must be so old (as member of a clean 
binary system) that gr has had the time to make many 
revolutions: 

5/3 / /Lf \ -2/3 

ta.e» — ^2xl0°yr| | ( — ) .(14) 



27r 



lO^days^ 

— The angular velocity of the periastron advance, ujr, 
has to be appreciably lager than the Galactic rotation 
(see Eq. (8)). 

In the latest pulsar catalog (Taylor et al. 1993, Taylor et 
al. 1995) we find more than 40 pulsars which are member 
of a binary system. Table 1 shows all small-eccentricity 
binary pulsars with P^ / e > 10^ days^. 



Table 1. List of test systems. The number given for the ec- 
centricity is the measured eccentricity plus the 2(7-error. The 
values for e were obtained by using nip = I.SMq and the upper 
limit for nic. 



Pulsar 


e < 


Pb 


P't/e> 


e > 






[days] 


[lO^days^] 




B0820-I-02 


0.012 


1232.5 


12.8 


6.9 


J1455-3330 ^ 


0.00017 


76.2 


3.4 


0.038 


J1640-I-2224 " 


0.00080 


175.5 


3.9 


0.21 


J1643-1224 ' 


0.00051 


147.0 


4.2 


0.25 


J1713-I-0747 


0.000075 


67.8 


6.1 


0.031 


B1800-27 


0.00051 


406.8 


31.9 


1.4 


B1953-I-29 


0.00033 


117.3 


4.2 


0.082 


J2019-I-2425 


0.00011 


76.5 


5.3 


0.037 


J2229-I-2643 "'^ 


0.00026 


93.0 


3.4 


0.053 



The binary pulsar PSR B0820-I-02 appears to have a 
hot, rather young white dwarf companion with an age of 
2 X 10^ yr (Koester et al. 1992). Thus Eq. (14) is not sat- 
isfied for PSR B0820-F02 for reasonable values of M. The 
binary pulsar PSR B1800— 27 shows the highest P^ /e and 
consequently would be at present the best test laboratory 
for a violation of the SEP. The period of the periastron 
advance is of the order of 3 x 10^ years. Unfortunately the 
evolutionary history of this pulsar is unclear. Like PSR 
B0820-I-02, which shows a similar position in the P-P di- 
agram (see Fig. 2), this binary system in its present form 
could be young, i.e. one cannot exclude that mass transfer 
took place within the past 10® years, which would lead to 
a violation of Eq. (14). This is supported by the fact that 
the characteristic age of the pulsar is only 3 x 10® years. 



^ see Lorimer et al. 1995 
" see Wolszcan 1995 
"'^ see Camilo 1995 



of the last mass transfer. Thus, for the rest of this paper, 
we shall handle this test system with some reservation. 

All the other pulsars in Tab. 1 are recycled millisecond 
binary pulsars with a white dwarf companion. All these 
systems are well away from the spin-up limit (see Fig. 2) 
and thus are old clean systems that clearly satisfy Eq. 
(14). 
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Fig. 2. P-P diagram for the pulsars of Tab. 1 



In Eq. (10) we need the total mass M and the orbital 
inclination i which is equivalent to the knowledge of 
and mc, see Eq. (13). The restrictions for these quantities 
are listed in Tab. 2. As mentioned before, in case of PSR 



Table 2. Limits on the parameters m^, mc, and i. For PSR 
J1713-I-0747 one gets further limits from the observation of the 
Shapiro delay caused by the companion 



Pulsar 


mp 


[Mq] 


mc 


[Mq] 


i 


[deg] 


J1455-3330 


1.0 


.. 1.8 


0.27 


.. 0.35 


0. 


. .90 


J1640-I-2224 


1.0 


.. 1.8 


0.31 


.. 0.40 


0. 


. .90 


J1643-1224 


1.0 


.. 1.8 


0.30 


.. 0.39 


0. 


. .90 


J1713-I-0747 


1.2 


.. 1.8 


0.27 


.. 0.34 


74 


. . .90 


B1800-27 


1.0 


.. 1.8 


0.36 


.. 0.47 


0. 


. .90 


B1953-I-29 


1.0 


.. 1.8 


0.29 


.. 0.37 


0. 


. .90 


J2019-I-2425 


1.0 


.. 1.8 


0.27 


.. 0.35 


0. 


. .90 


J2229-I-2643 


1.0 


.. 1.8 


0.28 


.. 0.36 


0. 


. .90 



J1713-I-0747 one was able to measure the Shapiro delay 
caused by the companion, which can be used to put further 
limits on i, mc, and (Camilo et al. 1994). 



Given the orbital period Ph, the eccentricity e, the total 
mass M , the orbital inclination i, and the Galactic accel- 
eration g, we can use Eq. (10) to determine the probability 
for a violation of the SEP, i.e. we calculate the area in the 
(^-Sl-parameter space where |A| is below a given value. We 
shall use the approximation ~ G in Eq. (10). Using 
Tab. 2 and assuming an error of 25% in the distance for 
the binary pulsars we obtain Tab. 3 (see also Wex 1994). 
Since all these systems represent statistically independent 



Table 3. Probability for the violation of the SEP being below 
a given |A| for eight test systems 



|A|< |A|< |A|< |A|< |A|< 



Pulsar 


1% 


0.5% 


0.4% 


0.3% 


0.2% 


J1455-3330 


0.82 


0.32 


0.00 


0.00 


0.00 


J1640-I-2224 


0.72 


0.37 


0.08 


0.00 


0.00 


J1643-1224 


0.73 


0.40 


0.23 


0.00 


0.00 


J1713-I-0747 


0.92 


0.82 


0.77 


0.61 


0.00 


B1953-I-29 


0.74 


0.44 


0.00 


0.00 


0.00 


J2019-I-2425 


0.80 


0.55 


0.43 


0.21 


0.00 


J2229-I-2643 


0.83 


0.24 


0.00 


0.00 


0.00 


El 


1.00 


0.99 


0.91 


0.69 


0.00 


B1800-27 


0.95 


0.92 


0.90 


0.87 


0.80 


S2 


1.00 


1.00 


0.99 


0.96 


0.80 



tests for the SEP we can give, besides an individual prob- 
ability for every binary pulsar, a total probability for the 
validity of the SEP denoted by S in Tab. 3. Since the 
evolutionary history of PSR B1800— 27 is uncertain (see 
previous section), we give two total probabilities, Si and 
S2, where the results obtained from PSR B1800— 27 are 
excluded in the case of Si and included in the case of S2. 

It is important to see that the restrictions on A do 
not contain any assumptions on the initial eccentricities 
of the binary systems. In this paper we made a worst case 
analysis, i.e. we assumed cr to be such that the largest 
possible is realised in the observed binary pulsars. If 
we would have a given probability distribution for cr we 
could get a better upper limit on A. 

In terms of a possible selection effect when choosing 
the test systems one has to assume that all the binary 
pulsars with a small /e were born like this. A fact that 
is very well explained by the eccentricity-period relation 
(Phinney 1992). This re lation explains the existence of 
an initial eccentricity which is different from zero for all 
the neutron-star white-dwarf binary pulsars and explains 
all the observed eccentricities for our test systems. 

If we express Eq. (1) in the field-theory based frame- 
work of Damour and Esposito-Farese (1995) and constrain 
the Nordtvedt parameter rj by Eqs. (2) then we find for 
the pulsar 

\ = hp[^ + C)+0{4). (15) 



(~ —ES™/mP) of the pulsar. For a medium equation 
of state one finds 

&p«1.03c2, Cp 0.21 nip/ Mq. (16) 

Since E^™ /rriiP ~ —10"'* for the white dwarf companion 
we have A ~ Ap, and thus Tab. 3 gives us limits on Ap. 
Using Tab. 3 and Eq. (15) one finds with 95% confidence 
for nip « IAMq: 

|e/2 + C|^5% (using Si), (17) 

|e/2 + C|^3% (using S2). (18) 

In Damour & Esposito-Farese (1995) one finds with 
90% confidence an upper limit of 1.5% for |Ap| . This 
result is solely based on an (obviously less conservative) 
analysis of the uncertain test system PSR B1800— 27 by 
Arzoumanian (see Damour & Esposito-Farese 1995). 

Bethe and Brown (Brown & Bethe 1994, Bethe & 
Brown 1995) give an upper limit for neutron-star masses 
of 1.56M0. This limit is supported by the assumption that 
a low-mass black hole was formed in SN 1987A. If we use 
this value as an upper limit for nip in Tab. 2, then we find 
with 95% confidence respectively 4% and 2.5% as upper 
limits for \e/2 + C\- 

5. Conclusions 

In this paper we presented new limits on the violation of 
the SEP in strong field regimes. The limits where achieved 
by analyzing long-orbital-period low-eccentricity binary 
pulsars. We made the conservative assumption that the 
mass of the pulsar lies between IMq and I.8M0. The 
mass of the companion was restricted by results of stellar- 
evolution scenarios of binary stars. For PSR J1713-I-0747 
we were able to use extra limits on the mass of the pul- 
sar and its companion that came from the observation of 
the Shapiro delay in the timing data. As a result of the 
increased number of suitable test systems we can exclude 
the violation of the SEP being above 0.5% with more than 
95% confidence. 

We were able to present limits on the tensor-multi- 
scalar theory of Damour and Esposito-Farese. For the two 
parameters of the (post) ^-Newtonian approximation of 
the tensor-multi-scalar theory we found the (conservative) 
limit |e/2 -|- C| ^ 5% with a confidence of 95%. Including 
the results of PSR B1800-27 lowered the limit to 3%. 

Using I.56M0 as upper limit for the masses of the 
pulsars, leads to a limit of 4%, respectively 2.5%, for \e/2+ 

C|. 
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